Microscopic mechanism for cold denaturation 
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We elucidate the mechanism of cold denaturation through constant-pressure simulations for a 
model of hydrophobic molecules in an explicit solvent. We find that the temperature dependence 
of the hydrophobic effect is the driving force/induces/facilitates cold denaturation. The physical 
mechanism underlying this phenomenon is identified as the destabilization of hydrophobic contact in 
favor of solvent separated configurations, the same mechanism seen in pressure induced denaturation. 
A phenomenological explanation proposed for the mechanism is suggested as being responsible for 
cold denaturation in real proteins. 
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Under physiological conditions, proteins adopt a 
unique three-dimensional (3D) structure It is max- 
imally stable at about 17°C and becomes unstable thus 
denaturing the protein at both high (~ 60° C) and low 
(~ -20°C) temperatures @, H, [|. The latter phe- 
nomenon is called cold denaturation, where the protein 
unfolds and thereby increases its entropy, which in turn 
is accompanied by a decrease in the entropy of the en- 
tire system. This counter-intuitive behavior has been 
experimentally verified [H, [f| but has remained a sub- 
ject of controversy @,|3l, since a satisfactory microscopic 
explanation for this phenomenon has not yet emerged. 
Resolving cold denaturation microscopically would facil- 
itate understanding the forces responsible for the struc- 
ture of proteins, and in particular the role of the complex 
hydrophobic effect 

In the case of diluted proteins, hydrophobicity is con- 
sidered as the main driving force for folding and unfold- 
ing Q . Consequently, different classes of models describ- 
ing hydrophobicity through varying explicit models have 
been used to study cold denaturation 0, H, H, [Io| ■ One 
such class 0, 0] considers the different energetic states of 
shell water, i.e. water molecules neighboring the protein, 
in a lattice. A more realistic water model [8] supports 
this view, as water-water hydrogen bonding among shell 
water has been found to increase at low temperatures and 
to correlate with cold denaturation. Meanwhile, another 
class of models suggests that the density fluctuations of 
water are responsible for cold denaturation [ll|, [l2| • De- 
spite the lack of consensus in the explanation of cold 
denaturation, the solvent is widely accepted as the key 
player. This is also supported by the fact that denatura- 
tion also takes place under pressure [HI, 0] • By focusing 
on the transfer of water molecules to the protein interior, 
pressure denaturation has been explained through the 



destabilization of hydrophobic_contacts in favor of sol- 
vent separated configurations 
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This destabilization 
has been verified using different water models [HI]. 

In the present work, we examine/unravel/reveal the 
microscopic physical mechanism behind cold denatura- 
tion. To this end, we consider the two-dimensional 
Mercedes-Benz (MB) model to describe water molecules 
in the solvent and a simple bead-spring model for the 
protein. The MB model reproduces many of the proper- 
ties of water |16[, including the temperature dependent 
behavior of the hydrophobic effect [l7j| . Our molecu- 
lar dynamics (MD) simulations of the MB model pro- 
vide a simple microscopic picture for cold denaturation 
in terms of changes in hydration: at low temperatures 
water molecules infiltrate the folded protein in order to 
passivate the "dangling" water-water hydrogen bonds 
(H-bonds) found in shell water. At the same time, hy- 
drophobic contacts are destabilized and an ordered layer 
of water molecules forms around the protein monomers 
such that they become separated by a layer of solvent in 
the cold denatured state. Hence, increasing pressure and 
decreasing temperature destabilize hydrophobic contacts 
in favor of similar solvent separated configurations. We 
expect that this aggravated destabilization of hydropho- 
bic contacts at high pressure explains why the transition 
temperature for cold denaturation increases with increas- 
ing pressure [2| . Here, we study cold denaturation at the 
equivalent of ambient pressure. 

As in water, the interaction between the MB molecules 
is given by a sum of hydrogen bonds and van der Waals 
bonds. The directionality of H-bonds is accounted for by 
three arms separated by an angle of 120°. This inter- 
action has maximal strength when arms of neighboring 
molecules are aligned. If is the distance vector be- 
tween the center of mass of molecules i and j, and r"i a 
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is the distance vector between the center of molecule i 
and the extremity of arm a, then the interaction energy 
is given by: 
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where and Rh are the binding energy and the equi- 
librium (reference) length of the bond, respectively. The 
constants <tr and <jg are attenuation parameters of the in- 
teraction. Equation (1) favors configurations where the 
distance between molecules i and j is Rh, one arm of 
molecule i is aligned with the line joining the two centers 
of mass, and the same for one arm of molecule j. The 
van der Waals interaction is described by a Lennard- 
Jones (LJ) potential V ww with binding energy e ww and 
equilibrium length R ww : 
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The LJ potentials are shifted so that the force becomes 
zero at the cut-off distance R c = 2.5Rh [Hj]. We use 
the parameter set which has been studied extensively by 
Silverstein et al. [l(|: eh = 1-0, Rh = 1-0, <jr — ag — 
0.085, e ww — 0.1, and R ww = 0.7. The total interaction 
energy Vij between two water molecules is given by the 
sum of Eqs. Q] and [2 

Here, we set M w = 1 for water. To mimic the distribu- 
tion of mass in water, 1/10 of the total mass of a water 
molecule is located at each arm's extremity and the ex- 
tremity of an arm is located at a distance i? a rm = 0.36Rh 
from the center of mass [l9j |. This defines the angular 
momentum of the water molecule. 

Energies, distances, and time are given in units of eh, 
Rh, and t = y/ (. W w I ^'lwR W w j respectively. To model 
the protein, we use a bead-spring model: monomers 
which are adjacent along the backbone of the pro- 
tein are connected to each other by springs, and non- 
adjacent monomers are connected by a shifted LJ po- 
tential. The LJ potential is described by a binding en- 
ergy e mm = 0.375 and distance R mm - The equilib- 
rium length and stiffness of the spring are i? S pring and 
-^spring = 2(456e mm /i?^ m ). This corresponds to twice 
the stiffness of the LJ potential. Monomers are set to 
be ten times heavier than water molecules. The interac- 
tion between monomers and water molecules is given by 
a shifted LJ potential with binding energy e wm — e ww 
and equilibrium length R wrn . 

When the side-chain of a hydrophobic amino acid 
is exposed to the solvent, the liquid surrounding the 



side-chain assumes a cage-like configuration 20] in or- 
der to minimize the amount of broken H-bonds of wa- 
ter molecules. This configuration has a low entropy and 
proteins minimize their free-energy by burying these hy- 
drophobic amino-acids in their interior. To reproduce 



this, wc choose R. 



spring 



= 2.0 and R„ 



= 0.9 such 



that monomers can be surrounded by a layer of water 
molecules when exposed to the solvent. To allow for the 
formation of a dry protein core, we use R mm — Rwm, 
though 2R mm > i? spr ing to avoid the backbone from in- 
tersecting itself. Taking these restrictions into account, 
we choose R mm = 1-1. 

Having defined the interaction between the different 
particles we now perform MD in the isothcrmal-isobaric 
ensemble. Constant pressure is achieved using the An- 
dersen extended method 2l|. To suppress oscillations 



of the simulation box, the canonical equations of motion 
are replaced by a Langevin stochastic process [13] imple- 
mented using the simplectic algorithm [23| . For the mass 
Q and the friction constant 7y of the piston acting on the 
simulation box, we use Q = 0.054/i?^ w and lv = 0.5. A 
parallelogram with equal sides and defined by an angle 
of 120° is used for the simulation box. This geometry re- 
tains the periodicity of a crystal made of water molecules 
through the boundaries. Periodic boundary conditions 
are implemented using the minimum image convention. 
For the Langevin equations describing the motion of par- 
ticles, we use the friction constant 7 _1 = 0.93r o . The 
noise term in the Langevin equations of motion is given 
by the fluctuation-dissipation theorem. Pressure is set to 
0.2 in units of thjR\- At this pressure, the MB model 
reproduces water-like anomalies seen at ambient pressure 



17j and hydrates non-polar molecules in a realistic man- 
ner [HI ■ The simulation box is packed with 512 molecules 
comprised of a 10-monomer long protein and 502 water 
molecules. To represent the solvent in its liquid state, we 
use temperatures ranging from 0.145 to 0.25 in units of 
e H . 

The system was initially equilibrated at a temperature 
of 0.25 for 5000 time steps, followed by a data collec- 
tion period of 50000 time steps. The temperature was 
then lowered and the equilibration-collection cycle was 
repeated. This cooling procedure was repeated until the 
lowest temperature was reached. Four samples with dif- 
ferent initial conditions were prepared using this proto- 
col and the distribution of the protein's radius of gyration 
Rg [HI was computed. To obtain equilibrium properties, 
the final configuration at each temperature was used to 
extend the simulation time until the distribution of Rq 
of the four samples converged within a root-mean-squarc 
value of 0.02. 

Let us now move on to describe the results. In Fig. [TJ 
we show the equilibrium distribution of Rq averaged over 
the four samples at three different temperatures. An ini- 
tial decrease in temperature, from 0.25 to 0.21, shifts 
the peak of the distribution to a lower value. Therefore, 
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FIG. 1: (Color Online) Normalized distribution of the radius 
of gyration Rg at three temperatures: T — 0.25, T = 0.21 
and T = 0.17. Inset: The temperature dependence of Rg of 
the protein. 

in hot water, proteins favor more compact configurations 
when the temperature of the system is lowered. However, 
a further decrease of temperature results in completely 
different behavior: as the temperature decreases from 
0.21 to 0.17, the peak shifts to a larger value indicat- 
ing that in cold water proteins become less compact for 
decreasing temperature. This behavior in hot and cold 
water is shown systematically in the inset of Figure [TJ 
which depicts the temperature dependence of Rq- The 
radius of gyration decreases as temperature decreases to- 
wards 0.21 - in hot water the protein folds as temper- 
ature decreases. Below that temperature, Rg increases 
monotonically as temperature decreases - in cold water 
the protein unfolds as temperature decreases. These two 
types of behavior are characteristic of heat and cold de- 
naturation of real proteins and are in line with previous 
studies 0,1113. 

The parabolic-like shape of Rg (see the inset of Fig. Q} 
cannot be mapped into a model with local monomer- 
monomer interactions only Q. To study the role of wa- 
ter, we show in Fig. [2] the average H-bond energy per 
water molecule for shell and bulk water. The energy of 
shell water averaged over the different configurations is 
higher than the energy of bulk water at high tempera- 
tures. This changes gradually as temperature decreases 
such that the creation of shell water becomes energeti- 
cally favorable at low temperatures. Therefore, when a 
protein is immersed in cold water it releases heat to form 
the shell, while in hot water it absorbs heat. These fea- 
tures are again characteristic of cold and heat denatura- 
tion of real proteins 0| • In the inset of Fig. [2] we show the 
energy absorbed by the system to create the shell around 
the protein. The absorbed energy is defined as the dif- 
ference in H-bond energy between shell and bulk water 
multiplied by the average number of molecules forming 
the shell. The absorbed energy decreases monotonically 



FIG. 2: (Color Online) Hydrogen bond energy per water 
molecule for shell and bulk water. Inset: Absorbed energy 
to accommodate the protein at different temperatures. The 
shell is defined by water molecules whose distance to the pro- 
tein is less than 2.5 in units of R^. 




FIG. 3: (Color Online) Characteristic configurations of a pro- 
tein in cold water (T = 0.15 and T = 0.17), at an intermediate 
temperature (T = 0.21), and in hot water (T = 0.25). The 
distance of highlighted (shell) water molecules to the protein 
is less than 2.5 in units of Rh- In cold water, the monomers 
are typically surrounded by clathrate-like cages. 



with decreasing temperature and becomes negative below 
some T indicating heat release. 

Characteristic configurations of the protein at different 
temperatures are shown in Fig. [3] In cold water (upper 
panels), the solvent forms a cage around each monomer 
of the protein, i.e. monomers are surrounded by an or- 
dered layer of water molecules. Molecules forming the 
cage are strongly H-bonded to each other and therefore 
have a low energy. At T = 0.21, the protein favors com- 
pact configurations. Water molecules close to the protein 
have at least one non-saturated H-bond which is pointing 
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towards the protein. When the temperature is increased 
to T = 0.25, most monomers are in contact with the sol- 
vent. The solvent forms either incomplete cages around 
monomers, i.e. cages which do not surround monomers 
from all sides, or they correspond to particles which are 
weakly bonded to the other solvent particles and are thus 
energetic. The crossover behavior of shell water shown in 
Fig.[2]is therefore characterized by the formation of cage- 
like configurations at low temperatures and the presence 
of dangling H-bonds at high temperatures. 

Configurations where monomers are separated by an 
ordered layer of solvent molecules have also been shown 
to become more stable, as temperature decreases, in 
models for the hydrophobic effect of methane- like solutes 
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26(| . Solvent layers around those monomer-pairs are 
highly ordered such that their formation decreases the 
entropy of the system 27|. Unfolding at low tempera- 



tures is therefore accompanied by a lowering in the en- 
tropy of the total system in accordance with experiments 
[3| , shell water molecules becoming more ordered as the 
protein becomes less ordered. This mechanism explains 
the counter-intuitive decrease in entropy during cold dc- 
naturation. The phenomenology is as follows. When 
non-polar solutes are transferred into water, the system 
relaxes by ordering those solvent molecules around the 
solute. This ordering has an entropic cost which is mini- 
mized by clustering non-polar solutes together, as this de- 
creases the amount of surface around solutes. As the tem- 
perature decreases below a particular value, the system 
rebuilds the ordered layer of solvent around non-polar so- 
lutes to saturate the dangling H-bonds left on the surface 
of clustered solutes - minimizing the enthalpy. Although 
hydrophobicity is not the only force responsible for the 
stability of proteins, the formation of a hydrophobic core 
plays the dominant role. 

In conclusion, we find that, at low temperatures, shell 
water forms hydrogen-bonds better than bulk water. Mi- 
croscopically this correlates with the presence of solvent- 
separated-configurations which accounts for the unfold- 
ing of the protein at low temperatures. The existence 
of such low energetic states for shell water at low T ex- 
plains why cold denaturation proceeds with heat release 
as opposed to heat absorption seen during heat denat- 
uration. Although here we studied cold denaturation 
in two dimensions, solvent-separated-configurations have 
also been shown to become more favorable as tempera- 
ture decreases in a 3D model for the hydrophobic effect 
[jjj ]. Therefore we expect that the results found in this 
work remain valid in 3D systems. Our results further sug- 
gest that cold and pressure denaturation could be studied 
under a single framework: a transition towards solvent- 
separated-configurations [l3| . 
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at T = 0.21. Therefore the entropy decreases when hy- separated configurations, 

drophobic contacts are destabilized in favor of solvent 



